ص۷۱۵ ام اتود ۵۲ افصتهول 


۹ 01/61۴ ۲۱088۵۵۵۵26۰ #ب-< 


۸ 1656210 
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۵ ۱0) ۵0۱۵ ۷ عظ) (0مه۵) 6۵ 0۲اهاصجصنی امصمت)هان مهم 66 مصمماع۱(۵۲ 
۵۸ ۵ 000۱۵۸7۵ 200 معتا۳۵ نو 0 
تتصوطاه ۱ ۲۱ مال0»‌زنم5 ,5 5 ,2172310000 1 11۰ 


صح ,نت120 ,ما1 ۵۶ نموه نصا رعصنمءصنعمط مصماوردمنظ ۵۶ احمنهمهن1 رعاهتا هتم ,۷۲56 -1 
صح 12012 ,متتطه1 ۵۶ اهتنا رعصتهمصنع مظ عصماووومزظ ۵۶ اصمعصصانهمع(1 ,۳۲۵]69501 -2 

صح ماه 1 ,مناج گم ونیته۲ نصا رعصتیمهء‌صتع صط عصصهاووعمزظ ۶ احممحصنهورهن1 ررمووع]۲0ظ ۸۵5[5]201 -3 
(2011710.26.11) ۵ 5660100 :1تمصظ عمط به عصت طمموعتمن -۶) 


۶۰ قفنطا عبات 0 ۲۱0۱۲ 


امصمتاهدمرجمی مط) عمتمم1۳۵۲۵1 .2022 .۱۱2192201 .۲۱ مج رتامال‌منلعگ .گ ٩9.‏ وبا رتتا2200 ۳2۲ 06-07-09 :156661۷60 
معط ممنصتاون0) مفصح موم لوط ۵ ممتانجم صمتاهاتاه۷ عمط امطصمت ما تماحانصنه 02-1-90 :16۷1860 
۳۵۲6120(۰ ع1) .131-146 :(2) 12 «تمصتطمه]۱۷ امتبلنه‌تتعه ۵۶ لمصتامآ صمتامنالمظ 07-12-09 :۸6660160 


101: 1022067/12 0 


(0 

6 ۱۵ 12010۲ م1 2 ور فصمتاتجمی صمصصصمتصه تقصماصا مطا عصتاامفجمع ,عععبامط ولمم «الهء‌تلتهعمو مه عامماوع11۷ ط[ 
عهع 40ص عسح‌هموصها نج مطا ملاناً معمتاتومی متتمطامومصاج مطا عماجم .عمتالهباعیه تتقطا امبمتی فصج ات تامنالمتم منهج 
معط ده فادها تقاصمصتهمظ .عممتاتجمع عصاتا مطا ۷۵محرصا صوع وعفامط تناو معلناً وععومره 0عوم1مصعتصهی صا ماج اصهعطهم 
بکاصتموم عماجم ۵۶ حمتهاتصنا ما مل مج ماماهم‌ع‌صهمالقطه میع وعفبامط وفلنامم هه مماوع۷ن1 مه فصمتاتجمی متمطام‌ومصطاه 
۵ 1۶ 1)1085جمی حصیتاصمجمه صرح متتحهم‌صصم) ستق مطا ۵۶ ممتاهااتگ .عیمه تهاصمصتهمنه همه فصمتاتصمی متقصتاه ماماه]وصتا 
۶ عصمتح‌اتهناً ما مصصمم‌تع۲ن 1۵ 05مطامصه (نات) ومتصمطرنا متا اهممتاه‌ نمی ها دممتتامومر فرصم له رالماتستاصتا 
لمصمتامعصظ ۶و معصهد اون مطه عصتام‌تهمنس 4صه فتورلفصه متتاعصه‌نهم ۵۶ ممتاتلاطاه اف۷2 عقط 0مطاعصط فتط1 .عادعا متصمصتهمرنه 
0 نا مصلوت 0ماهاتاصته ۱۷۷۵۲۵ عوبامط تواتاهس ج ۵۶ جمتباطاتتاعنل زاته‌ماع۷ 20 معا بنج مطا بجام‌نهعومر ونطة مر و .عتعامصه ت22 
۰ همه موجه مهتم مه ممتاهاتامه۷ عتة مطا 10 ممتات‌صمی ادها فطه ۵منطم2 


۷۵۵۵05 20 عامزه/۷]۵ 
اعد عط1 ,0مطعمصه مج مهبم ]رطاصصجن ممتونا 2060ع1ه م۷2 مقبامظ بوداوم گم 061مصط تقعتتامصصمعع فطا رتول‌باای مومت مط مر 
ط صصماطامتم عط 017۵و ما عمجم ۳/۵۲۵ عاصعصصماه 166550 ,عالتاعع1 عطا ما ع‌صتلتمععخ .۵0۵ص موله که ناک عمصعم1206۵۵ 

20 2606012016 ۰ 

عط 0تعصا تعلعصها صیاصم‌صصمص فطع متقاصصته ما 0مامعامو ۳2۵ ممتاهتاوع (5لنه۴) وماما م27 106-27072260معرم مط[ 
ص ممتام‌حصتاوعه صتقصط م1 .عصمتاه»ءنام‌مره توص عم واع0مهصط ممصماهمایجع 0عفجا اومص فط ۵۶ مجرم و1 ۵06ص ععع1 فظ1 ,عوبامظ وقلتا0ط 
عط صتطا حمتافنوع تعفصعته تقعط اصمتمصهتا ۸ .صوتصم ۱۱۵۳ و 0تسا مط فص فصه ماماتعوم»م‌صصممصا ود ملع مطع عقط ود [ممصط وتطا 
صمتامنلجهه مطا وعنامصا تقطا ودمعمزج کصمع‌جم060-مصصتا 2 وعماتجموع0 فقطا معتتههم‌جصها) تنج مطا امعم ما 0متام9۵16 7/28 صتقصمل با 
۵6 ۱۵۴۵۵۵۲۵۲۵ تتعطا 0صرع وعبامط 24 عصتست ر«الفامصتمون 0متاعهعصه ۱۷2۵ صمتاتصمی بونهل‌صتامها فطع ال .عصها ممتامع۲جمع 200 
0 ۱/۵۲۵ ۱۵0615 2200241621 م1 ,مق 0عمرم7۵1ع0 فص 6۵ 0متامرصره صرح قام0مصط اهمتاممصم‌طاممط تمممتور فطل عصتوهه 2۵06160 
نج ما عطنتاعومصظ رها راتم‌جمدنهون ۷۵102/604 وه تمتمانته 0عمم0۵7۵ مط1. ,عته۵ز ]۳1,7۳۲ - ۷۲5 ۸۸5 ما وادنامعصه)نصدزه 
۰(کاصلمح 13) وصمتاهع10 0م1لمعمو فصو ۵۶ عتنااجتهم‌صع] 


وی وزرا 0صه فالباجع1 

مج 60عاعجمج عطا معع۵)۲ه مبله۷ ۹۳۴ مطا کقط 5۵ رمهتاهعه امامح]فنامه ۵0مزصه 20۵061 مطا عقطا مهتاعصمصعل فااتاوع مظ[ 
عط داز متحهم‌جها سنح مطا امنلعتم ۱0نامی تمعاحصتی عط صرح 1.29 6۵ 0.405 ۵۶ ععصعد مطا صد م72 ممتتاجتهم‌هها تنج 2۲۵0060 
ومد بولامم ۵ عصتااه‌طجمی مصصتا- ۲۵2 مطا عم عماقانته 0ماهلنل۷2 ما مویا ما اصاتعووم و از رعع۵10تم1 .0621665 0.6 ۵۶ 260017207 
-16 صح )تاهج صرح معتعتهم‌جصعا عنة مطا ما واتعصعع۲0ع)عظ طعونط عطا رعتلتاعع1 مطا ما عصتلتمععظ .۳۵6896 ممتاهلتاصه۷ عط ه2تصصتامه ما 
له مطا بصمتات20 ص .عععبمط وفلنامج ۵۶ مممتلجعم1 مامتنه۷ 26 ممتاباطاتتاوتک متتحمم‌صها نع مص صد ۵0 موهای مه عمصمتعنل معتوع 
۵6 1 ,ععتامط متام ۵۶ واصتمم تحصصمی مطا صا «للفتهعووه رمفامط دمم ۵ عصهام اصمتعلنل معط له متا ام ۷۵۵ ]۷۵ 
۶ صمتاهاتاجه۷ عطا 0۷۵تحرصصا م) ۱960 ۱۲28 ۲ماهاتاصنه عط رمتمععط۲. ,عل۷2 060«مجصمرمعع۲ فطع صقط) معطع‌نط و قط ونعد 1 صقطا تمطونط 
عنه مطا عامعَ2 واصه۷ عاماصاً بنج عمط ۵۶ ماعصه عمط عقطا 0عاجعن4صا کیان 0مزه فممتاقامصته منامتیع ۵۶ مالناوع مط1- عففافظ تولتاوم فط 
جمده وومتمح م۷۵ عامه۷ ]فاص تج مه صمطاه صمنصها معمص منم مماتاماتمافتل ما۷ 24 معتمتمم‌مه] تنج مطا رفظ .عمصملناماننم 
نامه ۱.6 ۷۵متمحصز ما 0معصقده 2۵ وماهع عصتاعتیه معط ۵۶ ماعصه مطا فصه 0عصوومن ۷۷۵۲۵ وعاقع 6۷ مص5۵ ,۰]۳۵۲۵]01۲۵ ,تقطان 
له مط 0ج ۷۵0 مرح که۷ عقبامط توتاممر عط ۶م صمتاتصمی طمت)هاناصه۷ معط یصمتاهمتمصه طمناو برظ عقبامط نام مط ۵ صمتاتطمع 
6 086 ۲حصصتتن مص صد ۵0نموهان هط صقطا رمتصنا عمج متع۱7 وفبامط 260تصصتاورن معط ص ممتاماماتطاعنل متاجهم‌جعا 24 ما۷۵ 
6 رلهتامعمو۲۵ ,ی مر 0:46 فصح 0.23 مج (0* 24 ما 18) 16 297 م1 291 ۵۶ معممد عطا ها مهن« بواتهمله۷ فصح مسطهتمم‌ه) عند 
۰ ۵0۱ ۲0۶ ممتات‌جمی 1۲۵۵0060 عص کم مته عمیاله۷ 


«متعسلصم) 
عطا صا مامع۲۷ وومتمه ۵۶ عمتاهء‌تاوموه .صمتاناطاتطاعنل له عمط هم امعم اصهم‌تلتصونه 2 20ظ فلمع۷ فطا ۵۶ ملعصه عصتجهمه م1 
0 وه تماهاناصته 0۵۷۵10۵60 مط1 .عیااهمم‌جصع مه ما۷ تنج ۵۶ ممتاماماتطاعنل صمتصنا 0 160 عفبامط تام ام قللهس »0و 

۰ نا00 امتاتاقومی 4صه صوتععل م6 م2602 280 ۵610110۵066 


صمناحاتاصه ۷ ,ععجامط متام متا اه لممزمصل رومتصصهم ون فتباط تمممتاهات مهم :0۲۱۲۵۲8 > 


نشریه ماشین‌های کشاورزی 


7 ...مجح [//۱)05۰ظ 
که 


جلد ۰۱۲ شماره ۲. تابستان ۰۱۴۰۱ ص ۱۳۱-۱۴۶ 


توسعه شبیه‌ساز رایانه‌ای به‌منظور مدیریت بهتر شرایط تهویه‌ای مرغداری و بهینه‌سازی تولید 
حامد فرزانیور ؛ سیدصادق سیدلو هریس ا حبیبه نعلبندی؟ 


تاریخ دریافت: ۱۳۹۸/۰۴/۱۵ 
تاریخ پذیرش: ۱۳۹۸/۰۹/۱۶ 


چکیده 


تهویه یک متغیر کلیدی برای بهینه‌سازی شرایط محیط داخلی مرغداری است که با کنترل عواملی مانند درجه حرارت» سرعت هوا و غلظت گازهای 
موجوه می‌تواند از لفات طیور جاوگیزین کند. در حال خاضر نستورالممل‌های دقیقن در رسای طراحی و سباعت سبالن‌سای مرغتنازی و سیستم‌های 
تهویه با هدف کنترل دقیق شرایط محیطی توسعه نیافته است. در این تحقیق در یک مرغداری طولی و با تهویه مکانیکی فشار منفی توزیم دماو 
الگوهای جریان هوا با استفاده از دینامیک سیالات محاسباتی شبیه‌سازی شد و با استفاده از شبیه‌ساز توسعه‌یافته راهکارهایی برای بهبود سیستم تهویه 
ارایه شد. در توسعه شبیه‌ساز مذکور همه شرایط مرزی از طریق اندازه‌گیری‌های تجربی در طول یک دوره ۲۴ ساعته به‌دست آمد. اعتبارسنجی نتایج 
شبیه‌ساز» در ۱۳ نقطه آزمون شد و مشخصه‌های آماری بیان کننده دقت پیش‌بینی مانند 181۷97 بین دمای اندازه‌گیری شده و پیش‌بینی شده در 
محدوده ۰/۴۰۵ تا ۱/۲۹ قرار داشت و شبیه‌ساز قادر بود تغییرات دمای نقاط مختلف سالن مرغداری را با دقت ۰/۶ درجه پیش‌بینی نماید. نتایج شبیه‌ساز 
نشان داد که در مرغداری مورد مطالعه پروفیل دما در سطح استقرار طیورء حدود ۱۸ درجه در نقاط مختلف سالن تفاوت دمایی وجود دارد و سرعت هوا در 
نواحی مرکزی سالن بیش از یک متر بر ثانیه می‌باشد. بنابراین در مرحله بعدی از شبیه‌ساز توسعه‌یافته برای بهبود تهوبه سالن مرغداری مذکور استفاده و 
تغییراتی در محل استقرار ورودی های هوا و سرعت هوای اعمال و در نهایت با اصلاحات انجام شده توزیع دما در مرغداری نسبت به حالت اولیه 
یکنواخت‌تر شد و دمای هوا در محدوده ۲۹۱ تا ۲۹۷ کلوین (معادل ۱۸ تا ۲۴ درجه سانتی‌گراد) و در حد مطلوب قرار گرفت. هم‌چنین پس از اصلاحات 
انحام شده» سرعت هوا در اغلب نقاط سالن در محدوده‌ی ۳ : ۰/۴۶ متر بر ثانیه قرار داشت که در محدوده توصیه شده می‌باشد. بنابراین برای 


انجمن مهندسان مکالیک اپران 


طراحی و ساخت مرغداری‌ها و نصب سیستم‌های تهویه» شبیه‌ساز توسعه داده شده دارای عملکرد و دقت خوبی است. 


واژه‌های کلیدی: توزیع دم تهویهه سرعت هوا؛ شبیه‌سازی» مرغداری» یکنواختی 


مقدمه 


گوشت مرغ و تخم‌مرغ به‌عنوان منابع غنی از پروتئین» در طول 
سال‌ها به‌طور وسیعی در تغذیه انسان مورد توجه قرار گرفته است 
به‌طوری که در بعضی از کشورها که کمبود مراتع طبیعی مشهود است 
گوشت طیور به سرعت جانشین گوشت سایر مها شده و طی چند 
دهه گذشته مرغداری‌های طیور گوشتیء به یک صنعت بزرگ و 
سودآور تبدیل شده است. پیشرفت‌های اخیر در تولید گوشت مرغ و 
پرورش جوجه‌های گوشتی» حاصل اصلاحات ژنتیکی طیور استفاده از 


۱- دانش‌آموخته کارشناسی ارشد مهندسی مکانیک بیوسیستم» گروه مهندسی 

بیوسیستم. دانشکده کشاورزی, دانشگاه تبریزء ایران 

۲- استاد گروه مهندسی بیوسیستم. دانشکده کشاورزی, دانشگاه تبریزه ایران 

۳- استادیار گروه مهندسی بیوسیستم. دانشکده کشاورزی» دانشگاه تبریزء ایران 

(*+- نویسنده مسئول: 2010۴ ۵ ٩61160100‏ :اتقح۴) 
0 1022067/27 :10۲ 


جیره‌های غذایی کنسانتره» بهبود کنترل و پیشگیری از بیماری‌هاء 
آقدامات ایمنی و زیستی و هم‌چنین استفاده از فن‌آوری‌های نوین 
برای کنترل شرایط محیطی سالن‌های مرغداری می‌باشد که همه این 
افداماک عیب آفرایشی سیام و توش راید کته اس علی رش 
تمام اقدامات صورت گرفته تلفات حیوانی هنوز در مرغداری‌ها بسیار 
بالا است. در سیستم‌های تولید متراکم طیور تراکم بالای جوجه‌ها در 
سطح سالن» موجب کاهش جربان هوا و تجمع گازهای مضر در 
نزدیکی جوجه‌ها و کف مرغداری شده و نیز از میزان انتقال گرمای 
اضافی بدن طیور به محیط بیرون کاسته می‌شود که اين امر موجب 
بروز تنش گرمایی و مشکلات تنفسی و زیستی در جوجه‌ها می‌شود. 
در اثر فعالیت‌های زیستی طیور و نیز فعالیت باکتری‌هاء میزان اکسیژن 
در دسترس جوجه‌ها کاهش و درصد رطوبت. غلظت آمونیاک و 
گازهای مضر متصاعد شده از فضولات طیور افزایش یافته و کیفیست 
ای سالن مرغذاری اف پیدا می‌کند که با تاکافی بودن خریان هن 
در سطح طیور عوامل مذکور موجب افزایش تلفات حیوانی می‌گردند 


فرزانپور و همکاران. توسعه شبیه‌ساز رایانه‌ای به‌منظور مدیربت بهتر شرایط تهویه‌ای مرغداری... ۱۳۳ 


.]0 ۶ ۳۲6۵0065 :2001 .]۵ ۶ 2327 ۷ :1997 ,.۵1 6 5طمحصطصزه) 
(2006 بتااع ۵۱2۲0 0صه صفع۸1 :2003 


تونایی در به‌دست آوردن شرایط زیست‌محیطی مطلوب در 
سیستم‌های پرورش متراکم و مدرن حیوانات و به‌ویژه طیور بیشتر 
وابسته به نحوه طراحی و اجرای سیستم تهویه در اماکن نگهداری 
حیوانات است به‌طوری که با تهویه مناسب و بهینه می‌توان از میزان 
تلفات در سیستم‌های پرورش حیوانات و به‌خصوص طیور کاست. 
مشخصه‌های زیست‌محیطی فضای داخلی مرغداری‌ها و گاوداری‌ها از 
قبیل دماء رطوبت و غلظت گازهای اتمسفر سالن مرغداری» تابع 
الگوهای جریان هوا هستند. الگوهای جریان هوا پیوند اساسی بین 
محیط بیرون و محیط داخل ساختمان را برقرار می‌کنند. بنابراین برای 
پاسخگوپی به نیازهای برنامه‌های عملیاتی و مدیریتی در مرغداری‌ها 
درک درست از اصول جریان هوا به‌منظور حصول مقادیر مطلوب و 
توزیع مناسب آن ضروری است (2003 بصن 4صه 0). 

مطالعه تجربی الگوهای جریان هوا در سالن‌های مرغداری به 
دلیل وسعت سالن‌ها و تعدد دریچه‌های ورود و خروج هوا و نیز تاثیر 
متقابل عوامل روی یکدیگر, فرآیندی سخت. پیچیده. زمان‌بر و نیز 
هزینه‌بر است (1999 ,۲6۵06568 :2002 ,۷۷۵۱۱66۲ 20 فعاعقدان). 

به همین دلیل استفاده از تکنیک‌های جدید برای مطالعات مشابه 
مورد توجه قرار گرفته است. غلامرضایی و همکاران ( 00127622 
7 .4 6) از شبکه عصبی به‌عنوان یک ابزار مدل‌سازی 
غیرخطی برای پیش‌بینی دما و رطوبت در سالن مرغداری استفاده 
کردند. باغانی و آق‌خانی (2018 ,تعفط‌لطعه 4مه نعقطعد) با 
استفاده از مدل رگرسیونی نرخ انتشار گاز آمونیاک در سالن مرغخداری 
را با دقت بالا پیش‌بینی کردند. 

دینامیک سیالات محاسباتی (0۳۳) نیز در دهه‌های اخیر برای 
مطالعه شرایط تهویه‌ای مرغداری مورد استفاده قرار گرفته است. 
امروزه با استفاده از «۳ به‌کارگیری مدل‌های پیچیده‌ی ریاضی بر 
پایه معادلات دینامیک سیالات» روش‌های عددی و الگوریتم‌هایی که 
رفتار سیال را در فضاهای سه‌بعدی شبیه‌سازی می‌کنند. امکان‌پذیر 
شده است. در طول سالیان متمادی توسعه روش‌های مدرن محاسباتی 
همانند 0۳1 همراه با افزايش قدرت محاسباتی به‌عنوان نتیجه‌ای از 
توسعه فن‌آوری ساخت سخت‌افزارهاء استفاده از این روش‌ها را بسیار 
متنوع» دقیق و استانخ نموده است (2002 ,۹:۵2 0صد صعط)) 
به‌عنوان نمونه روش‌های عددی امکان آزمودن و ارزیابی ایده‌های 
طراحی را قیل از اجرا و یا ساخت مجموعه مرغداری فراهم می‌کند که 
یک مزیت فوق‌العاده محسوب می‌شود. 0۳1 از مشکلات و 
ناسا زگاری‌های مربوط به متشابه‌سازی اجباری عملکرد تهویه‌ای 
سالن‌ها به دلیل کاهش مقیاس» جلوگیری می‌کند زیرا که با استفاده 


هرن فتباظ تقممتامانام‌هون -1 


از د1) می‌توان در مقیاس کامل شبیه‌سازی را انجام داد. با استفاده 
از 27 می‌توان کنترل کاملی بر شرایط مرزی داششت و به راحشی و 
بسیار موثر تاثیر پارامترهای مختلف را مورد مطالعه قرار داد 
(2012 ,مع1ه10ظ 24 تصمحصه). تحقیقات زیادی با استفاده از 
0 در راستای پیش‌بینی الگوهای جریان هوا در مرغداری‌ها و 
دامداری‌ها انجام شده است که در اغلب این مطالعات شرایط آب و 
هوایی تحت شرایط یکنواخت شبیه‌سازی شده است. به‌طوری که لی و 
همکاران (2007 ,.0 2 6.]) در یک دامنه محاسباتی سه‌بعدی» 
پدیده‌های انتقال جرم و گرما را در یک مرغداری با تهویه طبیعی و در 
شرایط محیطی یکنواخت» شبیه‌سازی نموده و جریان هوا در داخل 
مرغداری را بررسی کردند. سئو و همکاران (2009 .6 2 860) با 
لحاظ دمای تولید شده توسط طیور در سیستم تهویه طبیعی» مدل 
توسعه‌یافته برای شبیه‌سازی شرایط محیطی مرغداری را ارتقا داده و 
به نتایج دقیق‌تری دست یافتند. در این مطالعه شرایط محیطی شامل 
دمای هوای محیط بیرون» دمای هوای داخل سالن مرغداری و 
گرمای تولیدی توسط طیور ثابت در نظر گرفته شد. تحت شرایط قابت 
دمایی» تحقیقات مشابهی نیز توسط ساير محققین انحام شده است 
(2012 ,.۵ 6۶ 2 :2012 ,اه ۶1 ۷]0۹]212(). 

روجانو و همکاران (2015 24 2014 ,.۵1 21 مصهژم؟) نیز طی 


تحقیقی با استفاده از (۳) محیط داخلی مرغداری مجهز به سیستم 
تهویه طبیعی را با در نظرگرفتن گرمای تولیدشده توسط طیوره گرمای 
حاصل از کود جامد و نیز گرماء رطوبت و دی‌اکسید کرین منتقل‌شده 
بین محیط داخل و بیرون مرغداری» با استفاده از مدل دو بمدی 
شبیه‌سازی کردند که نتایج منجر به افزايش صحت پیش‌بینی‌ها شد. 
مدل (۳۲) توسعه‌یافته در این مطالعه شامل پیش‌بینی حرارت 
رطوبت و دی‌اکسیدکرین تبادل شده بین محیط بیرون و داخل 
مرغداری بود. تمایز اصلی این مطالعه در نظر گرفتن گرمای منتشر 
شده توسط جوجه‌ها و کود جامد موجود در بستر مرغداری بود که 
به‌طور قابل‌توجهی پویایی شرایط محیطی را با توجه به تراکم بالای 
طیور تحت‌تاثیر قرار می‌دهد. با این وجود در مطالعه مذکور از متوسط 
دمای هوای ورودی در یک ساعت به سالن برای شبیه‌سازی استفاده 
شد. تحقیقات مشابه دیگری نیز توسط سایر محققین انحام شده است 
۱۷2۵060 20 ۷۷۵۲۱۵۷ :1994 ,اه 6 00۳۷۵۲۵۲ 2۵ ۷) 


:2005 ,اه 6 ۱۵0 :1997 ,اه 6 نامه :1995 
627 :2010 ,.]۵ ۶۶ 0۲۲۵۲ :2008 .اه ۶۶ 12865-۷1021 


(2014 .اه 61 12122 

نتایج حاصل از همه تحقیقات انجام شده نشان‌دهنده قدرت 
در شبیه‌سازی و پیش‌بینی شرایط داخل مرغداری‌ها می‌باشد 
که می‌توان از آن با دقت قابل قبولی برای اصلاح و مدیریت بهتر 
مرغداری‌ها به‌منظور کاهش تلفات حیوانی استفاده کرد. 

هدف تحقیق حاضر کاربرد دینامیک سیالات محاسباتی در 


۴ نشریه ماشین‌های کشاورزی. جلد ۱۲ شماره ۲. تابستان ۱۴۰۱ 


شبیه‌سازی نحوه توزیع دما و بررسی الگوهای جریان هوا در یک 
مرغداری می‌باشد تا با استفاده از تتایچ شبیه‌سازی بتوان پیشنهادهایی 
برای بهبود سیستم تهویه در راستای بالابردن راندمان تولید. از طریق 
اصلاح یکنواختی دمای محیط داخل مرغداری و به تبع آن یکنواختی 
شرایط محیطی داخل مرغداری ارایه کرد. در این شبیه‌سازی کلیه 
شرایط مرزی به‌صورت واقع‌گرايانه لحاظ شد به‌طوری‌که با 
اندازه‌گیری‌های میدانی و عملی دمای متغیر وافقعی سطوح مختلف 
شامل سقف. دیواره‌های جانبی» هوای ورودی از دریچه‌هاء دمای متغیر 
بستر مرغداری همراه با گرمای ساطع شده از سطح گرم‌کن‌های 
موجود در مرغداری در شرایط مرزی لحاظ شده است. 


مواد و روش‌ها 
مرغداری مورد مطالعه 

مرغداری مورد مطالعه از نوع طولی و با تهویه مکانیکی فشار 
منفی بوده و در موقعیت مکانی ۳۷/۹ درجه شمالی و ۴۷/۱۳ درجه 
شرقی قرار داشت. این نوع مرغداری» تیپ غالب مرغداری‌های چند 
استان ايران است که کنترل شرایط محیطی در آن‌ها نسبت به 
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مرغداری‌های عرضی نسبتاً آسان‌تر می‌باشد. هندسه و ابعاد سالن 
مرغداری در شکل 2-۱ نشان داده شده است. سالن مرغداری مذ کور 
به مساحت ۶۰۰ متر مربع در راستای شرق به غرب قرار داشته و 
دیوارها و سقف آن از بلوک‌های سیمانی ساخته شده بود. سیستم 
تهویه آن از نوع تهویه مکانیکی فشار منفی بود که برای ورود هوای 
تازه به داخل سالن مرغداری» روی هر کدام از دیوارهای شماره ۱ و 
۴ ۲ عدد دریچه با ابعاد ۴۰« ۱۰۰ سانتی‌متر و برای خروج هوای 
داخل سالن نیز ۴ عدد مکنده با قطر ۸۰ سانتی‌متر در دیوار شماره ۲ 
در ضلع غربی نصب شده بود. در زمان مطالعه» مرغداری دارای ۷۰۰۰ 
قطعه مرغ گوشتی ۴۰ و ۴۱ روزه و با جرم متوسط ۲/۵ کیلوگرم بود 
که به دلیل برودت هوا (آذر ماه) و برای انجام حداقل تهویه و حفظ 
دمای سالن مرغداری, فقط از سه دریجه در دیوار شماره ۱ و یک 
دریچه در دیوار شماره ۴ با زاویه بازشدگی ثابت ۱۵ درجه. برای ورود 
هوای تازه به داخل سالن مرغداری استفاده می‌شد (شکل 0-۱). تعداد 
دریچه‌های باز و محل آن‌ها و نیز مقدار بازشدگی آن‌ها صرفاً بر 
اساس تجربه مرغدار تنظیم می‌شد و هیچ داده‌برداری دمایی یا سرعت 
هوا برای انتخاب دریچه‌های باز و مقدار بازشدگی وجود نداشت. 


۷۲۲۵۱۱ 4 


۱۷۹3 


شکل ا- طرحواره مرغداری مورد مطالعه؛ 1 هیتر» 27 دریچه ورود هواه 1۹ مکنده» ۷۷ دیوا ۶ سقف 
۴ :5 :۷۲۷۵11 ۷۷۰ :صو۲ ومتامنا ۲۰ چاه ۷ ]علض تن ۷۰ :۲۱۵2/۵ ۲۱۰ :عوبامط توانامم 0منلبااو ۵ متاحتعتلدم ۲1۵۰1۰ 


هم‌چنین برای انتقال هوای داخل سالن مرغداری به بیرون نیز 
فقط ۲ مکنده (شکل ۱-) مورد استفاده قرار می‌گرفت. برای تامین 
گرمای مرغداری, از ۴ عدد گرم‌کن از نوع محفظه احتراق با ظرفیت 
حرارتی حدود ۱۰۰۰۰ کیل و کالری در ساعت) واقع در پایین 
دریجه‌های ورود هوا (شکل 2-۱ و ط) استفاده می‌شد. 
شبیه‌سازی شرایط محیطی مرغداری مورد مطالعه 
ایجاد مدل هندسی و المان‌بندی مرغداری 


(نسخه ۳-۲-۲) ترسیم و شبکه‌بندی دامنه محاسباتی در محیط آن 


انجام شد. ذر مدل مذکور صرفاً دریچه‌ها و مکنده‌های مورد استفاده 
توسط مرغدار در هندسه مورد نظر لحاظ و وضعیت موجود مرغداری 
ترسیم شد به‌طوری که فقط ۲ عدد از مکنده‌های انتهایی و نیز ۳ عدد 
از دریچه‌های ورود هوا (شکل 9-۱) توسط مرغدار مورد استفاده قرار 
می‌گرفت. جهت تسریع و کنترل بهتر محاسبات از شبکه‌بندی منظم 
استفاده شد. برای این منظور دامنه محاسباتی به دو قسمت محزا 
شامل فضای 2 و 9 تقسیم شد و برای شبکه‌بندی قسمت 2 از 
المان‌های مکعبی (۲1۵۵060721) و برای قسمت "از ترکیب 
المان‌های مکعبی و مثلشی (16172060721) استفاده شد (شکل ۲). 


دامنه محاسباتی اولیه دارای ۱۶۶۵۵۰ المان بود. 


فرزانپور و همکاران. توسعه شبیه‌ساز رایانه‌ای به‌منظور مدیریت بهتر شرایط تهویه‌ای مرغداری... 
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شکل ۲ المان‌بندی فضای محاسباتی مورد مطالعه 
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صتحصصمل 00و ۵۲ فترع آمصمتاهام۵ ۲۱۵۰2۰ 


معادلات حاکم 

پدیده‌های انتقال حاکم در محیط مرغداری شامل انتقال ممنتوم 
و انتقال گرما می‌باشند که برای شبیه‌سازی سالن مرغداری معادلات 
حاکم در محیط نرم‌افزار انسیس- فلوئنت با اعمال شرایط مرزی 
مناسب به‌صورت عددی حل شدند. 
انتقال ممنتوم 

معادلات حاکم در پدیده انتقال ممنتوم با فرض جریان غیر قابل 
تراکم و سیال نیوتنی به‌صورت روابط (۱) و (۲) بیان شدند که 
به ترتیب معادله بقای جرم یا پیوستگی و معادله بقای ممنتوم 
می‌باشند. 
)۱ ۲۰۷-0 


)+0 جج)+ ۳ ]۲۳-۲ )رم + ۷ ۳ 


استفاده از مدل آشفتگی در تحلیل مسئله 

جریان‌های تهویه مکانیکی به‌دلیل نرخ بالای جربان و نیز 
بهدلیل فعل و انعلات ناشی از نتقال گرما در میدان جریان, همراه 
با حرکت متلاطم هستند (2007 .۵1 2 0۲۵0). هم‌چنین در 
مطالعه انجام گرفته سرعت متوسط هوای ورودی از دریچه‌ها طی 
آزمایش‌های تجربی در حدود ۴/۲ متر بر ثانیه ثبت شد. بنابراین عدد 
رینولدز با استفاده از رابطه (۳) محاسبه شد و مقدار آن برابر با 


۶ به‌دست آمد که نشان‌دهنده متلاطم بودن رژيم جریان در 
محل دریچه‌های ورودی هوا است. در داخل سالن نیز به دلیل بزرگ 
بودن قطر هیدرولیکی مطمئنا رژیم جربان متلاطم پیش‌بینی 
می‌شود. بنابراین از مدل آشفتگی ع-* استاندارد که قابلیت زیادی در 
پیش‌بینی الگوهای جریان در جریان‌های تهویه مکانیکی دارد برای 
تحلیل مسئله استفاده شد (2009 ,اه ۶1 ظ۲01۲0). 
)۲ 

معادلات انتقال مذکور با استفاده از مدل آشفتگی ع-1 استاندارد 
به‌صورت روابط (۴) و (۵) بازنویسی شدند و رابطه (۶) به‌عنوان شرط 
اولیه برای تحلیل معادلات مورد استفاده قرار گرفت (-۸5۷5 
4 ,راظ۳۲[6). 


(۴) . قرم- ]۲ بتارم + 0 

ها ارم بطروهین +[۳۵ ۱۱/۵ +۲]0 نم 7 
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۶ شریه ماشین‌های کشاورزی. جلد ۱۲. شماره ۲. تابستان ۱۴۰۱ 


نمادها 
تا( 
2 سرعت زمان 
ره حص) بجازه0[ ۱۷ (5) عحطز[" 
دانسیته هوا : ضریب انتقال گرمای هوا 
۳5 (۳ 2 7اتعصول تتظ (1 ۳ ۷ 60000001 معط دنر 
8 عدد رینولدز ‌ ظرفیت گرمایی ویژه 
تحاحناه ۵/۵105 رای 1 «انمدمری احعط متلنهع5 
۷ جرم بدن طیور 5 منبع گرما 
(ع1) وعقصه وه تناو ( ۱0 
دما ویسکوزیته سیال 
7 (16) تمصع[ ِ راوخ ۵۲ ۳۲1010 
0 قطر هیدرولیکی مقطع سالن مرغداری 0 عدد ثابت برابر ۰/۰۹ 
(ه) عوامط پوتناوج ۵۶ تعامصحنل متلنیج۲۱۲۲ 9 0 601۲21 عاصع‌اعدمی ۷۲۵061[ 
۳ گرمای کل تولیدی توسط طیور در مرغداری ‌ عدد ثابت برابر ۱/۴۴ 
(۷۷) تاو عطا برد 0ع2تعصهع ۲۱۵۵۲ ۳ 44 ۵ ل2ناومع دهاجم ۷۲۵061[ 
1 دمای هوای سالن مرغداری ۰ عدد ثابت برابر ۱/۹۲ 
: (0) ععدامظ دمم 0۴ ماجومجه۱ تتظ 5 2 0 60۲21 عاحهاعدمی ۷۲۵01[ 
۳ انرژی جنبشی متلاطم ی عدد ثابت برابر ۱/۲ 
5 (ه که رععمجه معتاعصفا اصملناطنن۲ 5 3 0 60021 عاصماعد0ع ۷۲۵061[ 
17 ویسکوزیته متلاطم ی عدد ثابت برابر ۱ 
ی ۱۱۰۱ 1 0 69021 واصعاعدم ۷۲۵06[ 
انرژی اتلافی در اثر آشفتگی 82 ضریب تبیین 
(۶ که هه ماهنت ۵۶ مقر متلهوزووزظ طمتاحصتصماعل ۶ آموه‌تلوم) 
ی میانگین ریشه مربعات خطا 


1۹00]-113621-501012176 01 


برای حل معادلات مذکور شرایط مرزی زیر مورد استفاده قرار 

مرز ورود هوا: هوا از طریق دربچه‌های ورود هوا (1 ۰۷ ۰۷2 ۷3 
و ۷4) وارد مرغداری می‌شد که شرایط مرزی از نوع فشار ثابت (فشار 
اتمسفری) انتخاب شدند. زاوبه ورود هوا برابر میزان بازشدگی 
دریچه‌ها (۱۵ درجه) اعمال شد. 

مرز خروج هوا: مرزهای خروج هوا شامل مکنده‌های موجود در 
انتهای سالن مرغداری (۳1 و ۳2) بود که شرایط مرزی از نوع سرعت 
تعبین شد. متوسط سرعت هوای خروجی از مکنده‌ها طی آزمایش‌های 
تجربی به‌دست آمد. 

مرز دیوارها: سایر مرزهای موجود در دامنه محاسباتی شامل همه 
دیوارهای جانبی» سقف بستر و سطح گرم‌گن‌ها به‌صورت دیواره‌های 
بدون لغزش در نظر گرفته شدند. 
انتقال گرما 

تغییرات دمای هوای داخل سالن مرغداری منتج از گرمای ایجاد 
شده در اثر فعالیت‌های زیستی طیور گرمای ساطع شده از بستر 
مرغداری (حاصل از فضولات طیور)» دمای هوای ورودی به سالن و 
تغییرات دمای دیوارهای سالن در اثر تغییرات دمای محیط بیرون 


مرغداری و تابش آفتاب می‌باشد. بنابراین لازم بود معادله انتقال گرما 
در حالت ناپایا لحاظ و تحلیل شود. هم‌چنین پیش‌بینی می‌شد که در 
برجن :اضق سالن فرخدازی سرعت هو سر بابد و اسان کرت در 
این مناطق صرفا از طریق هدایت صورت گیرد. لذا فرآیند انتقال گرما 
به‌صورت ترکیبی از هدایت و جابه‌جایی درنظر گرفته شد (رابطه ۷). 
)۳ 7 + ۷۷۲ (0)00۷۲م + ۹ م0 

ی تا زور بای نع ادته فان سا استاده 
شدند: 

مرز ورود هوا: شرط مرزی در دریچه‌های ورود هوا (۸۷1 ۷2 
3 و ۷4) از نوع شرایط مرزی دمایی انتخاب شد. دمای هوای 
ورودی به سالن با آزمایش‌های تجربی در طول ۲۴ ساعت اندازه‌گیری 
و به‌صورت رابطه (۱۲) و تابعی از زمان مدل‌سازی شد. هم‌چنین دمای 
اولیه هوا برابر ۲۹۴ درجه کلوین به‌عنوان شرط اولیه لحاظ شد. 

مرز خروج هوا (مرز شماره ۲): از شرط مرزی شار همرفتی استفاده 

دیوار و سقف رو به آفتاب (1۳1۷4): برای دیوار و سقف رو به 
آفتاب مرغداری از شرط مرزی دمایی استفاده شد. دمای این مرزها 
نیز از طربق آزمایش‌های تجربی به‌صورت روابط (*) و (۱۰) 
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مدل‌سازی شدند. 

کف مرغداری: جهت اعمال گرمای ایجاد شده توسط طیور و 
بستر مرغداری از شرط مرزی شار گرمایی استفاده شد. با توجه به 
پراکنده بودن مرغ‌ها در کف سالن مرغداری و ارتفاع کم طیور» شار 
گرمایی در کف مرغداری در نظر گرفته شد. گرمای تولید شده توسط 
هر مرغ با استفاده از رابطه (۸) محاسبه شد (2015 ,.۵ ۶ م0طه[50). 
در این رابطه جرم بدن هر مرغ ۲/۵ کیلوگرم و متوسط دمای هوای 
محیط مرغداری ۲۱ درجه سانتی‌گراد در نظر گرفته شد. با در نظر 
گرفتن ۷۰۰۰ مرغ و مساحت کف سالن مرغداری» شار گرمایی برابر 
۵ وات بر متر مربع محاسبه گردید. گرمای تولید شده توسط 
فعل و انفعالات حاصل از بستر مرغداری نیز ۱۰ درصد گرمای حاصل 
از طیور در نظر گرفته شد. 
مه 

گرم کن‌ها (111 132 113 ۴14): برای اعمال گرمای تولید شده 
توسط گرم‌کن‌ها در فرایند شبیه‌سازی از متوسط دمای سطح 
گرم کن‌ها اسستفاده شند. گرم کن‌ها به‌ صورت دایم کار با فمای بت 
سطح برابر با ۲۳۹ درجه کلوین فرض شدند. جنس گرم‌کن‌ها از فولاد 
با دانسیته ۸۰۳۰ کیلوگرم بر متر مکعب. گرمای ویژه‌ی ۵۰۲/۴۸ ژول 
بر کیلوگرم درجه کلوین و با ضریب هدایت حرارتی ۱۶/۲۷ وات بر 
مترسدرجه کلوین لحاظ شد. 

سایر مرزهای دامنه محاسباتی: برای سایر مرزها (۷1 7۷2 
3 2؟) نیز از شرط مرزی دما استفاده شد. تغییرات دمای سطوح 
مذکور نیز همانند دمای سقف و دیوار رو به آفتاب مدل‌سازی شد 
(رابطه ۱۱). 


۸۵1 - 0۹ + 


حل عددی 

برای حل عددی معادلات ناوبر استوکس از الگوریتم حل مبتنی 
بر فشار استفاده شد. جهت تحلیل و حل مساله با استفاده از کوپلینگ 
فشار- سرعت استفاده شد. در این مطالعه. از روش فشار ضمنی پیزو 
(۲ماحهع0 ۵۶ عصتانامک طانه اهنا عسوعع۳۳) برای ارتباط 


1,0,0 


900000 


بین سرعت و فشار با توجه به ماهیت مسئله و حل آن در حالت ناپایه 
استفاده شد (2014 ,/۸(55-۳[۵60), 

شبیه‌سازی مرغداری طی ۲۴ ساعت کارکرد آن انجام گرفت. حل 
عددی معادلات با گام‌های زمانی ۱۰ ثانیه‌ای و در محموع 
شامل۸۶۴۰ کام با شناد حتاکقر. ۱۰۰۰ تکرار سرا هر گام زسانی 
راهاندازی و تنظیم شد. 

شرایط همگرایی: باقی‌مانده‌های مورد قبول برای حل معادله 
پیوستگی برابر ۰/۰۰۱ برای معادلات سرعت. انرژی جنبشی 
توربولانس و نیز نرخ تلفات انرژی جنبشی توربولانس برابر ۰/۰۰۰۱ و 
در نهایت برای معادله انرژی برابر ۰/۰۰۰۰۰۱ در نظر گرفته شد. 


جهت بررسی استقلال شبکه چهار نوع شبکه‌بندی به‌ترتیب با 
۸ ۰۲۶۶۵۸۵۰ ۶۲۱۳۱۲ و ۱۰۳۰۰۱۶ المان با شرایط مرزی 
یکسان انتخاب و تحلیل شده و نتایج دما استخراج گردید. بر اساس 
نتایج به‌دست آمده (شکل ۴) اختلاف بین سطوح سوم و چهارم کمتر 
از ۳ درصد بود. لذا برای کاهش ححم محاسبات 9 صرفه‌جویی در 
زمان تحلیل از تعداد ۶۲۱۳۱۲ المان برای تحلیل‌های نهایی استفاده 


آزمایش‌های تجربی 

دمای داخل سالن مرغداری همواره تابعی از دمای محیط بیرون 
مرغداری است که تابعی از فصل و نیز زمان روز متفاوت خواهد بود. 
دمای محیط بیرون مرغداری» گرمای منتقل شده از سقف و دیوارهای 
سالن و نیز دمای هوای ورودی به سالن را تحت تاثیر قرار می‌دهد. لذا 
برای افزايش دقت شبیه‌سازی لازم بود از مقادیر واقعی در شرایط 
مرزی مدل توسعه‌یافته استفاده شود. برای این منظور آزمایش‌های 
تجربی گسترده‌ای در ۲۴ام آذرماه سال ۱۳۹۶ انجام و تغییرات دمای 
مرزها در طول یک دوره ۲۴ ساعته. اندازه‌گیری و مدل‌سازی شد. 
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معطععجط گم تهطاججیت از 


شکل ۳- نتایج بررسی استقلال شبکه 
۷ 060606060 هد وعصه ۵ عناوم .3 ,۲12 


۸ شریه ماشین‌های کشاورزی. جلد ۱۲. شماره ۲. تابستان ۱۴۰۱ 


لازم به توضیح است شرایط مرزی از جمله دیوارها و سقف و نیز 
دمای هوای ورودی به سالن تابع زمان روز و تاریخ داده‌برداری خواهد 
بود. به دلیل متفیر بودن این داده‌ها در طول سال و برای امکان 
استفاده از شبیه‌ساز توسعه‌داده شده هميیشه لازم است که متوسط 
دمای دیوارها در یک بازه زمانی در ماه يا فصل مورد نظر و نیز 
متوسط دمای هوای ورودی و متوسط تغییرات دمای هوا در طول روز 
در همان بازه زمانی اندازه‌گیری و در مدل اعمال شود تا وضعیت 
پروفیل دمایی و سرعتی در سالن مشخص گردد. 

دمای دیوار و سقف رو به آفتاب: برای اندازه‌گیری دمای سطح 
دیوار و سقف رو به آفتاب (۱۷4 11) روی دیوارهای مورد نظر یک 
عدد ترموکوپل نوع ‏ (با رنج عملکرد دمایی بین ۲۰۰- الی ۱۳۵۰ + 
درجه سانتی‌گراد و حساسیت /۴۱۳۸۷) نصب و در مدت ۲۴ ساعت 
تغییرات دمای آن‌ها هر ۵ دقیقه یکبار به‌وسیله دیتالاگر (12ظ) ثبت 
شد. داده‌برداری از ساعت ۱۲:۳۰ بعد از ظهر شروع شد. داده‌های ثبت 
شده در محیط نرم‌افزار متلب به‌صورت تابعی از زمان مدل‌سازی شده 
و در شرایط مرزی اعمال شد. این نوع معادلات بر اساس منطقه 
جغرافیایی و زمان فصل باید تبیین و تعیین شود. 

دمای دیواره‌ها و سقف رو به شمال (۷71 7/2 ۷3 032: 
تغییرات دمای سطح دیوار ۷۷1 به‌عنوان نماینده‌ی همه دیوارها و 
سقف قرار گرفته در سایه نیز به روش مشابه در قسمت قبل. 
اندازه‌گیری و در نهایت تغییرات آن‌ها تابعی از زمان مدل‌سازی شد. 

دمای هوای ورودی به سالن مرغداری (۷1 ۰۷2 ۷3 و ۷4): 
دمای هوای ورودی به سالن مرغداری نیز به روش مشابه اندازه‌گیری 
دمای دیوارها و در دریچه‌های ورود هوا به داخل سالن (۴ دریچه در 
شکل 0-۱) اندازه‌گیری و مدل‌سازی شد. 

دمای سطح گرم کن‌ها (11 ۰112 113 ۳14): برای تعیین متوسط 
دمای سطح گرم کن‌ها از دو ترموکوپل استفاده شد. یک ترموکوپل 
روی سطح یکی از گرم کن‌ها نصب و دمای سطح آن در طول ۲ دوره 
روشن و خاموش شدن مشعل اندازه‌گیری شد. ترموکوپل دوم در 
مجاورت گرم کن نصب شد تا تعداد دفعات و مدت زمان روشن و 
خاموش شدن گرم‌کن در طول ۲۴ ساعت» مشخص گردد. در نهایت 
با توجه به تعداد دفعات روشن خاموش شدن گرم‌کن و هم‌چنین 
میانگین دمای سطح گرم‌کن در طول ۲ دوره روشن و خاموش شدن» 
متوسط دمای سطح گرم‌کن محاسبه شده و در مدل‌ها به‌عنوان منبع 
حرارتی اعمال شد که این امر برای ساده‌سازی مدل توسعه یافته بود 
در غیر این صورت بایستی تغییرات دمایی منبع حرارتی به‌صورت ناپایا 
تابعی از زمان در مدل‌ها اعمال می‌شد. 

سرعت هوای خروجی از مکنده‌ها (۳1 و ۳2): این سرعت 
به‌وسیله یک دستگاه سرعت‌سنج از نوع سیمداغ (1۳5-1341) 


اندازه گیری شد و سرعت هوا در ۸ نقطه در مقطع هر کدام از مکنده‌ها 
ثبت و میانگین آن‌ها به‌عنوان متوسط سرعت هوای خروجی در نظر 
گرفته شد که در طول آزمایش‌ها هميشه ثابت بود. 


اعتبارسنجی مدل توسعه‌یافته 

هم‌زمان با اندازه‌گیری‌های مربوط به شرایط مرزی داده‌برداری 
دمایی در نقاط مختلف سالن مرغداری و در شرایط تهویه واقعی برای 
اعتبارسنجی مدل توسعه‌یافته انجام شد. به دلیل بزرگ بودن سالن 
مرغداری, ساختمان به دو بخش ۸ و ظ تقسیم شد (شکل ۴) و دمای 
هوا در ۱۳ نقطه مختلف در طول ۲۴ ساعت به‌وسیله ترموکوپل 
اندازه‌گیری و ثبت شد. نقاط تعیین شده برای ثبت دمادر ۲ ارتفاع 
مختلف سالن شامل ارتفاع ۰۴۰ ۱۴۰ و ۲۰۰ سانتی‌متری از بستر سالن 
مرغداری قرار داشتند. برای اعتبارسنجی مدل توسعه‌یافته دمای 
پیش‌بینی شده توسط شبیه‌ساز در هر نقطه با متوسط دمای 
اندازه‌گیری شده مقایسه و با استفاده از شاخص‌آماری ۳ دقت 
مدل توسعه‌یافته بررسی شد. 

به غیر از نقاط نزدیک دریچه‌های ورودی و خروجی هوا؛ هميشه 
مقادیر سرعت هوا در مقاطع مختلف سالن بسیار بطنی و ناچیز 
می‌باشد که این مقادیر برای تخلیه گازهای سمی و غیرسمی حاصل 
از تتفس و فضولات کف سالن کفایت می‌کند فقط بایستی توزیع 
سرعت هوا در نقاط مختلف یکنواخت باشد. اندازه‌گیری سرعت‌های 
هوای کم با استفاده از ابزارهای موجود تقریباً دشوار و در مواردی 
دارای دقت بسیار پایینی است که در این تحقیق نیز چنین مشکلی 
حادث شد. بنابراین با توجه به این که پروفیل‌های سرعت و دمای هوا 
به همدیگر وابسته می‌باشند پس اعتبارسنجی مدل صرفاً بر پایه 
اندازه‌گیری‌های دمای هوا صورت گرفت. 


بهینه‌سازی مرغداری 

در پایان بر اساس نتایج حاصل از شبیه‌سازی مرغداری با۴ 
دریچه باز برای ورود هوا و ۲ مکنده در حال کار وضعیت مرغداری از 
لحاظ یکنواختی توزیع دما و نیز توزیع سرعت هوا در نقاط مختلف 
مرغداری مورد مطالعه قرار گرفت و برای ایجاد یکنواختی بیشتر در 
ففیت خذویه داغلن:(دما و سترغت هوا اصلاحای :از فلز اقرایفن 
تعداد دریچه‌های باز و مقدار زاویه بازشدگی آن‌هاء توزیع نحوه باز 
شدن دریجه‌های ۲۴ گانه و سرعت هوای خروجی از مکنده‌ها اعمال 
گردید به‌طوری که طرح نهایی اصلاح شده در بخش نتایج آورده شده 


است. 
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سانتی‌متری و نقطه ۱۲ در ارتفاع ۲۰۰ سانتی‌متری از کف مرغداری) 
حصم من 40۵ 0۲ اطعتمط :9 6۵ 1 واصتمم) ۲ماقلناصتی 0عمم7۵ع0 عمتامل‌نله۷ ما عاصتمم عمتاعععصه مسعهم‌صها 0 1/0620 ۲۱۵,4۰ 
(0صیامع مطا حمظ جصن 40 ۵۶ )طونم :12 عصام زلصنامع مط) مج ۵۶40 )طونم :13 هه 11 ,10 فاصتوم ج 0صنامتع فطا 


نتایج اندازه‌گیری شرایط مرزی 

برای افزايش دقت مدل توسعه داده شده. همه شرایط مرزی 
ب‌صورت تجربی اندازه‌گیری شدند. 

دمای دیوار و سقف رو به آفتاب (۷۷4 1): تابش مستقیم آفتاب 
به دیوار رو به آفتاب (۷4) و سقف رو به آفتاب (4181 سبب گرم 
شدن دیوار شده که اختلاف دمای آن با هوای داخل سالن مرغداری 


0) 


۳ 
تغییرات دمای سقف رو به سایه و سایر دیوارها به غیر از دیوار رو 

به آفتاب (۰۷/1 ۷2 ۷3 1۳2): اختلاف دمای دیوارهای قرار گرفته 
در طرف سایه با دمای هوای داخل سالن مرغداری» سبب شار گرمای 
منفی و هدر رفت گرما از این دیوارها به طرف بیرون می‌شد. تغییرات 
دما در سطح داخلی آن‌ها به‌صورت تابعی از زمان روز و به‌صورت 
رابطه (۱۱) مدل‌سازی و در شرایط مرزی در حل معادلات استفاده شد. 
(« 1+2 896*1077 5+ 7*۷1 7.003<10- 2 رون[ 
05 - ۷۲۹۳۴( ,0.7085 2< *18 


تغییرات دمای محیط بیرون: دمای هوای داخل سالن مرغداری 
به‌طور قابل‌توجهی متاثر از دمای هوای وارد شده به مرغداری بود. 
تغییرات دمای هوای ورودی به سالن مرغداری تابعی از زمان و 
ب‌صورت رابطه (۱۳) مدل‌سازی شد. 
1+0 *<0.0003668- 7*1 10*< 4.042 < ورن[ 


۱۲ 
1* - 0.8078, (۷۲۳ 2 6 


موجب ایجاد شار گرمای مثبت و یا منفی به سالن مرغداری شود. 
تغییرات دما در مرز ۷۷4 به‌صورت تابعی از زمان روز و در طول ۲۴ 
ساعت به‌صورت رابطه )٩(‏ و تغییرات دمای مرز 1*1 به‌صورت تابعی از 
زمان روز و به‌صورت رابطه (۱۰) مدل‌سازی شدند. شاخص‌های آماری 
و 18۷15۳ که بیانگر اعتبار مدل می‌باشند برای هر کدام از 
مدل‌های شرایط مرزی در زیر آورده شده‌اند. مقادیر آن‌ها بیانگر دقت 
بالای مدل‌های استقرایی توسعه داده شده می‌باشد. 
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(*<0.00006431) و 0.5157 + (0 *<0.00006431) ومع 0.09696 -292.3 < رم[ 
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نتایج اندازه‌گیری سرعت هوای مکنده‌هاو نیز دمای سطح 
گرم‌کن‌ها نشان داد که مقادیر آن‌ها به‌ترتیب برابر ۷/۲ متر بر ثانیه و 
۳ کلوین (معادل ۲۳۹ درجه سانتی‌گراد) بودند. 
اعتبارسنجی مدل توسعه‌یافته 

همان‌طوری که قبلاً ذکر شد اعتبارسنجی مدل توسعه‌یافته در ۱۳ 
نقطه انجام شد. مقادیر 1۷18۳ بین دمای اندازه‌گیری شده و 
پیش‌بینی شده توسط مدل در محدوده ۰/۴۰۵ ت۱۲ ۱/۲۹ قرار داشت 
(جدول ۱). مقادیر به‌دست آمده در این تحقیق بهتر از نتایج ارایه شده 
توسط روجانو و همکاران (2015 .۵1 ۶ 0حهزم۳) بود که در آن 
مقدار ]1۷ در محدوده ۰/٩‏ الی ۱/۹ گزارش شده بود. نمودار 
مربوط به دمای ۲ نقطه به‌عنوان نماینده‌ای از نقاط اندازه‌گیری شده 
در شکل ۵ آورده شده است که نشان می‌دهد تطابق خوبی بین 
داده‌های تجربی و شبیه‌سازی شده وجود دارد به‌طوری که در نقطه ۳ 
حداکثر حداقل و متوسط اختلاف بین دمای آندازه‌گیری شده و 
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شبیه‌سازی شده به‌ترتیب برابر ۱/۶۴ ۰/۰۰۱ و ۰/۴۹۷ درجه بود. در 
نقاط ٩‏ و ۱۰ نیز متوسط اختلاف دما به‌ترتیب برابر ۰/۶ و ۰/۷۵ درجه 
به‌دست آمد. تفاوت‌های مشاهده شده بین داده‌های تجربی و شبیه- 
سازی شده به دلیل محدودیت و صحت پیش‌بینی معادلات حاکم 
استفاده شده در بیان دینامیک انتقال جرم و انرژی توسط حیوانات بود. 
به‌عنوان مثال رطوبت به‌عنوان عنصری موثر در شرایط محیط داخلی 
مرغداری می‌تواند در مقدار تولید گرما موثر باشد ولی در معادلات 
حاکم برای شبیه‌سازی در نظر گرفته نشده است. هم‌چنین در فرآیند 
شییه‌سازی دمای گرم‌کن قابت در نظر گرفته شند ولی در شرایط واقمی 
گرم‌کن در دوره‌های زمانی معینی روشن و خاموش می‌شد. با این 
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وجود مدل توسعه‌یافته به خوبی قادر بود تغییرات دمای نقاط مختلف 
سالن مرغداری را با دقت ۰/۶ درجه پیش‌بینی کند و می‌توان از آن 
برای مطالعه شرایط محیطی مرغداری و مدیریت بهتر آن استفاده 
کرد. 
نتایج حاصل از شبیه‌سازی 

پس از اعتبارسنجی مدل توسعه‌یافته. شرایط تهوبه‌ای مرغداری 
شامل پروفیل‌های توزیع دما و سرعت هوا در نقاط مختلف سالن مورد 
مطالعه قرار گرفت که نتایج آن در ادامه آورده شده است. 
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شکل ۵- منحنی تغییرات دمای پیش‌بینی شده و اندازه گیری شده در نقاط مختلف سالن مرغداری در پرابر زمان 
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مصصا که ماع 


جدول - مقادیر 1۷1 بین دمای اندازه‌گیری شده و پیش‌بینی شده توسط مدل در نقاط مختلف سالن مرغداری 
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پیش‌بینی پروفیل سرعت هوای سالن مرغداری 

همان‌طوری که قبلا اشاره شد برای تهوبه مرغداری در فصل سرد 
از ۴ دریچه (شکل 0-۱) استفاده می‌شد که انتخاب تعداده محل و 
میزان بازشدگی آن‌ها بر اساس تجربه مرغدار صورت می‌گرفت و 
هیچ‌گونه اطلاعاتی در مورد یکنواختی توزیع سرعت هوا و دمای هوا 
در سالن مرغداری در دست نبود. نتایج حاصل از شبیه‌سازی شرایط 
محیطی مرغداری نشان داد که غیریکن واختی قابل‌توجهی از لحاظ 
سرعت هوا در نقاط مختلف سالن مرغداری وجود دارد به‌طوری که در 
برخی نقاط سالن سرعت هوا نزدیک صفر و در برخی مناطق بیش از 
مقادیر توصیه شده بود. در مقطع عرضی واقع در مجاورت دریچه‌های 
۱و ۴ (شکل 2-۶ و 0) بیشترین سرعت هوا در نواحی مرکزی 
مرغداری حدود ۲/۵ متر بر ثانیه مشاهده شد. در مقطع عرضی در 
مجاورت دریچه شماره ۳ (شکل 4-۷) نیز توزیع سرعت هوا یکنواخت 
نبود و در ناحیه 2-2 شکل 2-۷ سرعت هوا حدود صفر بود. 

توزیع دما و تجمع گازها در سطح فرارگیری طیور (ارتفاع ۰/۴ 


متری از کف سالن مرغداری» شکل ۸- ۵) از دیدگاه زیستی و 
سلامتی مرغ‌ها حائز اهمیت می‌باشد. حداقل مقدار سرعت هوا در 
شد که عدم تهویه مناسب در این نواحی موجب اجتماع طیور در 
مناطق مرکزی مرغداری شده و تهویه مرکز سالن را با مشکل مواجه 
قسمت‌هایی از نواحی مرکزی سالن و نواحی مجاور گرم‌کن شماره ۱ 
(۳1) مشاهده شد. سرعت بالای هوای سرد در سطح طیور موجب 
کاهش دمای بدن آن‌ها و مشکلات زیستی می‌شود. به همین دلیل 
سرعت هوای مجاز در کف مرغداری و در سطح طیور در فصول سرد 
و در زمان تهویه حداقل» کمتر از ۰/۲۵ متر بر ثانیه توصیه شده است 
(2007 ,.۵1 21 8ع). ملاحظه می‌شود الگوی قرارگیری دریچه‌ها در 
دیوره‌های جانبی بیش‌ترین تأثیر را روی یکنواختی جریان هو دارد 
که در مرغداری مورد مطالعه به درستی انتخاب نشدند. 
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شکل - بردار و پروفیل سرعت هوا (2 و 0) و پروفیل دمای هوا () در راستای عرض سالن مرغداری و در مجاورت دریچه‌های شماره ۱ و ۴ پس 
از ۱۰ ساعت 
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شکل ۷- پروفیل بردارهای سرعت هوا (ه) و دمای هوا (0) در مقطع عرضی سالن مرغداری و در مجاورت دریچه‌ی شماره ۳ پس از ۱۰ ساعت 
۲ 099۵ ناو ۵۴ ممتاعع5 ققمتم 2 (0) ماجهم‌ها عنع ۵۴ ماتلمتی هه (ه) واله‌ماع۷ تنج ۵۲ ملتلمیج فصح ما60 ۷ ۲۱۵۰7۰ 
فسامط 10 تمه 3 واصه۷ مط 


البته باید توجه داشت که باز کردن تعداد بیشتری از دریچه‌های 
ورود هوا به دلیل افزايش دبی هوای ورودی به سالن مرغداری موجب 
یکنواختی بیشتر سرعت هوا در سالن شده و با تهویه مناسب‌تر» موجب 
کاهش تلفات طیور می‌گردد که این مسئله در فصول گرم سال 
مشکل خاصی ایجاد نمی‌کند و صرفاً لازم است طراحی محل 
دریچه‌هاء زاویه بازشدگی و میزان دبی هوای ورودی به مرغداری به 
شکل صحیحی انتخاب شود. ولی در فصول سرد جهت حصول 


شکل ۸- پروفیل سرعت () و دمای هوا () در راستای طولی مرغداری در ارتفاع ۴۰ سانتی‌متری از کف سالن (سطح طیور) پس 


یکنواختی بیشتر سرعت هوا در مرغداری» استفاده از دبی‌های بالای 
هوای ورودی هرچند موجب کاهش تلفات طیور رشد بهتر و افزایش 
وزن طیور می‌شود ولی هم‌زمان مصرف انرژی بابت گرمایش سالن 
افزایش خواهد یافت. بنابراین ضمن این‌که طراحی, اجرا و تنظیم 
دریجه‌ها و مکنده‌های مرغداری باید بر اساس نتایج شبیه‌سازی باشد 


علاوه بر آن بایستی اثرات اقتصادی آن نیز مدنظر قرار گیرد تا هزینه 
تولید واحد محصول بهینه گردد. 
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فرزانپور و همکاران. توسعه شبیه‌ساز رایانه‌ای به‌منظور مدیربت بهتر شرایط تهویه‌ای مرغداری... ۱۳۳ 


پروفیل دمای هوای مرغداری 

پروفیل توزیع دمای پیش‌بینی شده توسط شبیه‌ساز در مقطح 
عرضی سالن مرغداری پس از سپری شدن ۱۰ ساعت از زمان شروع 
مطالعه (ساعت ۲۲:۳۰) در شکل ۰-۶ و ۷- 9 نشان داده شده است. 
ملاحظه می‌شود دمای هوا در مجاورت بستر سالن مرغداری به علت 
مجاورت با گرم‌کن و هم‌چنین مجاورت با محل استقرار طیور به‌عنوان 
منابع تولید گرماءبالاترین مقدار (هترتیب ۵۱۰ و ۳۱۰ کلوین معادل 
۷ و ۲۳۷ درجه سانتی‌گراد) و در مجاورت سقف مرغداری پایین‌ترین 
مقدار را دارد. لازم به ذکر است مقادیر دماهای گزارش شده وابسته به 
زمان و دمای هوای بیرون سالن بوده و با تغییر آن‌ها تغییر نمودند. 
بررسی‌های بیشتر پروفیل دمای مرغداری در راستای طولی نشان 
می‌دهد که دما در نواحی مرکزی و قسمت انتهایی سالن مرغداری 
(محاورت دیوار شماره ۲) که در نزدیکی مکنده‌ها قرار دارد نسبت به 
نواحی دیگر بالاتر است که علت آن نوع تهویه مرغداری و حرکت هوا 
از سمت دیواره‌های جانبی به سمت مرکز سالن و در نهایت حرکت آن 
به سمت انتهای مرغداری و به سمت مکنده‌ها می‌باشد. در طول این 
مسیرء تباال گرمای همرفتی بین هوله طیور, بستر و گرم‌کن‌ها موجب 
فزایش تدریجی دمای هوا در طول سالن مرغداری می‌شود. در مقاطع 
عرضی سالن بین دمای هوا در نواحی کناره دیوارهای جانبی و 
مرکزی سالن حدود ۱۰ درجه اختلاف دما وجود دارد. در ابتدا و انتهای 
سالن مرغداری نیز بیش از۱۴ درجه اختلاف دما مشاهده شد که مقدار 
آن دز ارفاخ‌هاین شاف ارت برد تایه شانییی فد تاسط سای 


محققین گزارش شده است .۵ 6 522722 :2012 ,.۵1 61 ۷05۲212) 
(2012 ,.۵1 ۶۶ عاصحصطعاعداظ ,2013 

توزیع یکنواخت دما در سطح طیور بسیار با اهمیت‌تر از سایر 
دهد. مطالعه پروفیل دما در سطح طیور نشان می‌دهد که 
غیریکنواختی قبل‌توجهی در دمای هوا در سطح سالن وجود دارد 
(شکل 0-۸). به‌طور متوسط ۱۸ درجه اختلاف بین دمای هوا در نقاط 
گرم کن شماره ۳ اختلاف دما به بیش از ۴۰ درجه رسیده است. هرچند 
غیریکنواختی توزیع طیور در سالن مرغداری در شبیه‌سازی انجام شده 
در نظر گرفته نشده است ولی تهویه ناکافی و بالا بودن دمای هوا در 
برخی نواحی سالن مرغداری موجب تجمع طیور در نواحی با شرایط 
مناسب‌تر از لحاظ دمایی مانند مرکز سالن خواهد شد که خود این امر 
موجب افزایش دمای این نواحی و غیریکنواختی بیشتر دمای هوا 
می‌گردد. 


اصلاح پروفیل توزیع سرعت و دمای هوای مرغداری 

همان طوری که نتایج بخشر قبل نشان داد توزیع سرعت و دمای 
هوا در تمام نقاط سالن مرغداری یکنواخت نبود. بنابراین تغییراتی 
جهت بهبود شرایط موجود روی مرغداری انجام شد (جدول ۲) که به 
نتایج هر مرحله در ادامه اشاره می‌شود. 


جدول ۲- تغییرات انجام شده در سالن مرغداری به‌منظور بهینه‌سازی توزیع سرعت هوا و دما 
مها 0صه چازهم(۷۵ تنج ۵۲ جمتاجاطاتاعتک مطا 0۷۵تصرصص مه عقبامط تناو قطا ۵۶ صمتاهع۱۷۲۵0۲11 -2 2016 ]1 


تغیبرات انجام شده 
6 0۱218266 
تغییر زاویه بازشدگی دریچه‌ها 


نتیجه به‌دست آمده 
تنااننت)۱ 


زاویه دریجه‌های ۰۷1 ۰۷2 ۷3 و ۷4 از ۱۵ به ۱۰ درجه کاهش یافت. 


صاصع۷ داعلصا تنج عط) ۶ ماعحضظ 
۱ 


باز کردن دو دریچه جدید 


قجع۷ ۵ ۱۳۵ عصت0ظ 


باز کردن دو دریچه دیگر و تغییر در زاویه 
بازشدگی دریچه‌ها 
۶ ۷2۱۵100 200 مامه۷ 070 عصتل خر 
کاصع۷ داعلصا نع عطا ۵۶ ماعص2 
کاهش سرعت در مرز خروجی (طرح نهایی) 
عطا صا واته‌ماه۷ ته مطاً هصتفدعت۵ع1 
رصعتععل لفصتل) )ما0 


تغییر در زاویه بازشدگی دریچه‌ها 
با توجه به این که در مجاورت دریچه‌های شماره ۱ و ۴ سرعت 


۰ 10 ۲0 15 1۳0 060۳62560 ۷۵۲۵ عاجع۷ ۷4 220 ۷3 ,۷2 ,۷1 0۶ لعج مط1 


با حفظ تغییرات مرحله ۸ دو دربچه دیگر (۷5 و ۷6) با زاویه بازشدگی ۱۰ درجه به سالن مرغداری افزوده شدند (شکل ۴). 
تا عطا ] 20060 7۷۵۲۵ (۷۵6 20 ۷5) داهع۷ ۱۷۷۵ ,1 ععهاه 0۴ «مناهتنه۷ عط) عصتصتهاصتقصط برظ 
۰( ۳1۵۰) 162۳66 10 ۵۶ ع1عصه عصتصهون برها 0256 
دو دریچه‌ی دیگر (۷7 و ۷8) با زاویه بازشدگی ۱۰ درجه اضافه شدند (شکل ۴) و زاویه بازشدگی دریچه‌های ۲ ۳ ۵ و ۶ از 
۰ درجه به ۱۵ درجه تغییر یافت. 

۰) 068766 10 ۶ه عاعصه ع‌صنصهمه ود فقتامظ پوتتاوج فص ۵ 20060 مهب (۷8 24 ۷7) داوه۷ 1۷۷0۵ 
6 15 ۲0 10 108 1۳60162560 ۱۷۵۲۵ ۷6 200 ۷5 ,۷3 ,۷2 0۶ ع1عصه عصنصهوه عط) 0ص (4 
با حفظ تغییرات مرحله ۲ سرعت هوا در مکنده‌ها از ۷/۳ متر بر ثانیه به ۶/٩‏ متر بر ثانیه کاهش یافت. 

0 08073 102660 ۷/۵8۵ 1۵0/6۲5 0۶ تزانه۷(0 نج ما ,3 معها؟ 0۶ جمتاهننه۷ عط) عصتصتهاصتقصط برظ 
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هوا در سطح طیور بیش‌تر از حد مجاز برای فصل سرد بود و به تبع 
آن دمای هوا در اين ناحیه نیز پایین‌تر بوده جهت بهبود این شرایط 


۱۳۳ نشریه ماشین‌های کشاورزی. جلد ۰۱۲ شماره ۲. تابستان ۱۴۰۱ 


زاویه بازشدگی دریچه‌ها از ۱۵ به ۱۰ درجه کاهش داده شد. این تغییر 
موجب یکنواختی توزیع سرعت هوا در سالن و یکنواختی توزیع دما در 
نواحی مجاور گرم‌کن‌های ۱ و ۴ در سطح طیور شد. با این وجود در 
نواحی مجاور دریچه‌های ۲ و ۳ و در سطح طیور تغییرات اعمال شده 
سبب افزایش دمای هوا در ناحیه مرکزی مرغداری نسبت به حالت 
اولیه شد. 
باز نمودن دو دریچه جدید (۷5 و ۷6) 

نتایج مطالعات نشان داد که وجود دو دریجه مقابل هم به دلیل 
ایجاد تلاطم, موجب اختلاط بهتر هوا در نواحی مرکزی سالن 
می‌گردد. بنابراین دو دریچه جدید (۷5 و ۷6) با زاوبه بازشدگی ۱۰ 
درجه در دیوار شماره ۴ و در مقابل دریچه‌های ۲ و ۲ ایجاد شد. با 
تغییرات انجام شده توزیع دما و سرعت هوا در مجاورت گرم‌کن‌های 
مربوط به دریچه‌های شماره ۲ و ۲ بهتر و ناحیه کور در مرکز 
مرغداری نسبت به حالت اولیه محدودتر شد. ولی در محاورت دیوار 
شماره ۳ و در مرکز مرغداری هنوز مناطقی با دمای هوای بالا و 
تهویه نامطلوب وجود داشت. 


باز نمودن دو دریچه‌ی دیگر و تغییر مجدد زاویه بازشدگی 
دریچه‌ها 

با توجه به مطلوب نبودن نتایج مرحله قبل, دو دربچه دیگر (۷7 و 
8 در دیوارهای شماره ۱ و ۴ با زاویه بازشدگی ۱۰ درجه باز شده و 
زاویه بازشدگی دریجه‌های ۲ ۲ ۵ و ۶ از ۱۰ درجه به ۱۵ درجه 
افزایش یافت. با تغییرات انجام شده یکنواختی توزیع دما و سرعت هوا 
در نواحی مرکزی مرغداری نسبت به مرحله قبل و نیز حالت اولیه 
مرغداری بهتر شد. ولی هم‌چنان در سطح طیور شرایط تهوبه دمای 
هوا مطلوب نبود. 

کاهش سرعت هوای خروجی از مکنده‌ها 


برای رساندن سرعت هوا به حد مطلوب و با حفظ تغییرات انجام 
شده در مرحله قبل» سرعت هوای خروجی از مکنده‌ها از ۷/۳ به ۶/۹ 
متر بر ثانیه کاهش داده شد. بر اساس نتایج حاصل.ء دمای هوا در 
سطح طیور افزايش يافته و به حد مطلوب رسید به‌طوری که در اغلب 
نواحی مرغداری دمای هوا در محدوده‌ی ۲۹۲ تا ۲۹۷ کلوین (معادل 
٩‏ تا ۲۴ درجه سانتی‌گراد) و سرعت هوا نیز در محدوده‌ی ۰/۲۵ تا 
۷ متر بر ثانیه قرار داشت. 


طرح نهایی 

با اعمال تغییرات مذکور طرح نهایی مرغداری به شرح زیر 

دریچه‌های ۱ ۴ ۷ و۸ با زاویه بازشدگی برابر ۱۰ درجه؛ 
دریجه‌های ۰۲ ۲ ۵ و ۶ با زاویه بازشدگی برابر ۱۵ درجه و سرعت 
مکنده‌های خروجی برابر ۶/۹ متر بر ثانیه (زشکل ). 

با اعمال تغییرات فوق, نتایج حاصل از شبیه‌سازی شرایط تهوبه 
مرغداری اصللاح شده نشان داد که میزان یکنواختی توزیع دما نسبت 
حالت اولیه مرغداری بیش‌تر شده و در دمای هوا در محدوده ۲۹۱ تا 
۷ کلوین (معادل ۱۸ تا ۲۴ درجه سانتی‌گراد) و در حد مطلوب قرار 
درد (شکل 4۱۰ پروفیل سرعت هوا در سعلح طیور نیز نان داد که 
قاط با سرعت هوای صفر و نقاط کور در مرغداری اصلاح شده نسبت 
به حالت اولیه کمتر است و هم‌چنین یکنواختی و دامنه اختلاف 
سرعت هوا در حالت اصلاح شده کمتر می‌باشد. در مرغداری اصلاح 
شده سرعت هوا در اغلب نقاط مرغداری در محدوده‌ی ۰/۲۳ متر بر 
انیه تا ۰/۴۶ متر بر ثانیه قرار دارد که به مقادیر توصیه‌شده توسط لی 
و همکاران (2007 ,.۸ 2 6ع]) برای سرعت هوا در طول فصل سرد 
و تهویه حداقل نزدیک است. 


شکل -٩‏ طرح نهایی سالن مرغداری 
۵ 00۱010 6ط) 0۲ صعتوعل لهصز۲ ۲1۵,۵۰ 


فرزانپور و همکاران. توسعه شبیه‌ساز رایانه‌ای به‌منظور مدیریت بهتر شرایط تهویه‌ای مرغداری... 
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شکل ۱۰- پروفیل دمای هوا در راستای طولی مرغداری اصلاح شده در ارتفاع ۴۰ سانتی‌متری از کف سالن (سطح طیور) پس از ۱۰ ساعت 
تمه 0صنامع عمط و200۷ مد 40 ,ععتامط نام 01160مجظ عط) ۵۶ ممتاععو لمو7۵عصها مطا که متنجوم‌جصم) بنه ۵۶ ۳۲۵1116 ۲۱۵۰10۵۰ 


تام 10 


نتیجه گیری در داخل مرغداری مورد استفاده قرار گرفت. نتایج حاصل از 


مه ۰ 


شبیه‌سازی‌های مختلف در جهت بهبود تهویه سیستم نشان داد که 


تهوبه یک متغیر کلیدی برای بهینه‌سازی شرایط محیط داخلی در 
مرغداری‌ها است و به‌طور مستقیم نرخ تلفات و بهره‌وری تولید را 
تحت تأثیر قرار می‌دهد. در حال حاضر دستور العمل‌های دقیقی در 
راستای طراحی و ساخت سالن‌های مرغداری و سیستم‌های تهویه 
توسعه نیافته و اکثر یبن پیشنهادهای ارائه شده براساس تحربه عمومی 
است. در این تحقیق برا الگوها بان هوا و توزیع دما 
ای ِ بررسی آلگو ِِ كِِ ِِِ ‌ ن مرغداری و در سطح طیور حاصل شد. 
یک مرغداری تیپ موجود در منطقه و به‌منظور بهبود یکنواختی توزیع 
دما و یکنواختی سرعت هوا مدل ۳۲) توسعه یافت. مدل توسعه‌یافته 
پس از اعتبارسنجی جهت پیش‌بینی الگوهای جریان هوا و توزیع دما 


زاویه بازشدگی دریجه‌ها بیشترین تاثیر را در اختلاط هوا دارد. 
هم‌چنین استفاده از دریچه‌های روبه‌روی هم در دیوارهای جانبی 
موجب یکنواختی بهتر توزیع سرعت و دمای هوا می‌گردد. بر اساس 
نتایج حاصل از شبیه‌سازی عددی. سالن مرغداری مورد مطالعه 
بهینه‌سازی و حداکثر یکنواختی توزیع سرعت و دمای هوا در سالن 


ز 3 ۱۰4 


۲5 ۷ 5- ۳۱۵۵8۲۰ 2014. 1۳601۷7 ۰ 

نام راتمدوری اصمعالنل صا 62760 وتملزمعه هم وتوراقصه عمجم مایت .2006 .باع0عقص‌امک .۱ هه ریق مطععاتاضظ 
393-0۰ :(4) 5 معصم1ع افص ۵ لقصرستاول صهنله)] .وععبامط 

۶ ۳۲۲6۵۲5 معط هه وعوبام ملظ وم ممتووتحطط فهین متعمصصصض ۵۶ تاو .2018 .تصمططي۸ ۲۰ ۱۷۲ 4صه ,.۷۲ ,تصه‌طعه‌ظ 
:۱01 ۳۵۳۹1۵۲(۰ ط1) .377-388 2(۰) 8 مدمه لحعتانه‌نبهد ۵۶ تمصدمل ۲ چم عم همه بانتصا! موم[ 
۰ 3 1-۱-۱( 

0 وعتصمطرول مت لمممتامانامهمی که ممتاهع۵11مظ .2008 10۳۲۵9۰ .0۰ بخ 24 رطام‌فحاظ .9 ,هتتهزنتان .ظ .۷ ,102 ۷-وعصع1ظ 
0 م۲۱6۵ فصماوروماظ .مصتلانه بقانمم لمتهتعصمم 0ع)ه1نام۷۵ اله)‌نصمطهعصظ ه ص «مائتته که ممتاعتلعتم فط 
105-۰ :(1) 

۲ .2012 .۰101۲6۵8 .0 .۸۵ 200 رتعلمانمدم۳۱ ۸ بطمعهلظ .9 ,ر0عمانا-فاهتهن ٩.‏ ۲۰ ,متتهزتنان .ظ ر.ظط رصاطهصصهافباظ 
:(5) 12 (ا286ظ) وع0قصعو .عصصع؟ بتلنامم صا فصمت)تصمی عتاقصصتاه 1۳000۶ ووعوفه م6 فصممطم‌عتا۱۵28ظ 1801610۵0۲۵1 0۲ حطع]5(٩‏ 
.-5752 

صمطمصت)هل رفصمتاناهه ما نامع 2 :فصماماميی اصمصصممت مه پتالناوظ .2002 .تملل2 ۷۷ ۷۷۰ هر هه ریک .1 روعاعقط) 
۰ ۷۵۲516۷ ]۲ 

لل ]موجه تمملصا ۵ مصتامو10 24 مصمتاهلنله۷ بممتاه۷۵ 10۲ 0۲06۵60۲۵ ۸ .2002 9۲۵0116۰ .۲ 200 ریما مصفطل) 
201-۰ :(2) 8 16862100 ۸۵ ۲۲۷ .2021۷565 

:۵ ۱01۲۳0 20 مطمتاهلنامه تنج رماع ممتاقتامه ۲۷ :2003 1۵۷۵۲۰ 1۱۰ 20 ,ما2 ۱۷۲ راعناطتمصصنهدر؟ .۲ ورگ ,۲6006 
:(3) 12 طمتهعومک و۲ 0عن1ومخ ه لمصتاهل م1 .عمصمصصموم بعلزميه مج فتاتابالع ۵۶ معصعلمصا فطا وم فا۵ع۳۲۶ 
.326-4 

ممتامصاوصمن) رعتمعظ مج واتتصتا ,تمصع .2017 تصطامعا هل بکا مه رتتع20 بک بلط بله2عتصدامدان 
-546 :(2) 7 «مصنطمه اه‌هد ۵۶ تفصیامل معط لمعبواه لمته‌تایه عملهونا القلظ وتنام مطا صا عصتافهدعع۳۲۵۲ 
۰ 1022067/12 :1۵۲ ۲۵۲۹1۵۴(۰ 1۶) .557 

12 46 1880106060124 .2014 .دقرم 06۵۵0۲۵ .۲ مه رحعطممهه ع6م بخ ۷۲ ب805۵۵0 مقعماماهی نک ریگ م2قکن 00۵24162 
۵۸ 0۱2 ,612ص متتحوم اع0 مععمته مه فهع4کناوون0 ععطامو. مقامگم20منمصعو. مباقفوظ . أع. ۱۷۲۵0111620168 
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۱۳۶ نشریه ماشین‌های کشاورزی. جلد ۱۲. شماره ۲. تابستان ۱۴۰۱ 


۲۵۵۵۱۸۵0 ۵۵10821 48 ۸۵۲0۴0۳۱24 67 2( 5۰ 

تمه ه که مناد ممتاهاتامه۷ مطا 1۵۲ پمجتتهعع بات ۵۶ ممتاح‌نله۳۲۷ .2007 .عصناگ .ظ .9 4صه ,موه .و ر.ظ بآ ,عع1 
.53-4 :(1) 41 مره :ناهام طمتهموم آمزبه‌لنمتتهم مهمرجز ,ععبامط تماتمتها 0ماقهاتامه۷ 

5۲ 00 2 وه (۲۲)) ومنمصمه(۱ ت۴۱ آمممتامنام‌جمن :1997 0۲ظ .0 هه رفصمططمع۷۷2 .۷ رعط0ل .12 .1 ریق مفتام۱۷۲1501 
نموه ه تمصدم وهصهاعطها .عمصتهانهط نموه ما ماقحصتاهمتمنصظر رممهصا مه ممتاهاتامه۷ ۵۶ وتورلقمه مطا 
۵0612006(۰) م51 

101 .1۰ .۳۲ 200 رع10ظ .ظ . رعظه/۲۱۳ .9 ,۲ ,ق۲۱۵۲ ۷۷۰ .9 ,960 با بل و60 .9 یک ,9018 .۲ .9 بععا .ظ ب1و.ظ مقلهاد۷]0 
۲۱۱۵0 عصلللنیاه عملزم:ه تفص ممتاطاتتونل مجهمجها عنم گم ممتاه‌اته ومتصمصول تس لمممتاهنام‌همن :2012 
.293-۰ :(4) 112 عط1عممم1عظظ فصصماووماظ ,ورد صمتاهاتامه۷ 0۲6۲ 

۱۵۵1۲ 0۴ ووممه۷ 611600 ممتتهاتامه۷ عظ1 عملووعوو۸ .2009 ماو ۷۷۰ 12۰ 200 رحجملله۲ یک راهن . .1 بطمات0 
فصماوءماظ .ومتصهصرل . مت تمممتامان موم . عصویا.. فصمتتصمی . 0مامصتصم . مص . بعلصتا.. فعصتل‌آتداط . ع1مم90ع11۷ 
78-۰ :(1) 103 عط۲1ع۳۴۴186 

امه ۵ ممتاه‌وتلدمی مملاهاتامع۷ مطا عصتفتتامن .2010 میاه ۷۷۰ .نا مه رعماله یک رطقعت) .لگ و.گ بطماع۱0 
4(۰) 45 )ممحصوم1ظ مصح عصتلاتباظ .تمصع محومط لمممجصهم هه تممص 40ع7معمرصصا عم فعصتل‌لنناه عل0م)9ع۷ 1 0ع)ه1ا۷۵ 
983-5۰ 

ومتصقول مب تمصمتمن موم ۵ مممتاهء‌نان۸۱ .2007 .۱2۵۵0 ۷۰ 220 رصملله۲ نک راصهتت رل تاک ۷۷۰ .م1 وگ بط0۲)0 
16001 ۲با0وعع10ظ ,۲۵۷1۵7۷ ۵ :4ص تحتانه‌تنوه معط مد فمصموروه ممتماتاصوب ۶ صعنوع مه عصذاع0مصه فطه مد (نات) 
.2386-14 :(12) 98 

۵و .2013 .۷۵1۵9006۵5 0 ۰[ ۳۲۰ هه ,ملمطمه]۱۷ 5۱۳۷۵ .۱۱ رقم 0۵11۷۵۲۵ .و یک رفصتاته فع۸۵0 ,۷ ,ین م۸ رگ ر2هته 
۰۵ 204 682201۷70 ۷1۵ وموبامط تملژمعه طا 0تون همم ماقااصته ما ومتصصمصول متا لمممتلمنا مهم 0 
159-۰ :(1) 16 0۵0011102 ۱۱۷۱۵۱۱۵۵06 مک 020ناهنتامه هن انا م)۷15عک .عصم)درد ممتاهاتامه۷ 9۵۵ امصصت 
مدمه ی ممتامج ۲۶۵0 تقصتصه 1۲ ماه ۲۷ رعصتیهعصلعصظ امعتانه‌تبوه که ۳۱2۵0000 010۴ .1999 .5 بطه‌وتع0عظ 
0 2 .2.3 ]27ظ رولقصتصض ۵۲ اممصطمتز صظ 2 عمامرم راصمصمنتاصظ مصح عصتویامز۲ >۲۸۷۵۵)۵0 1 )۲2۲ رقطلبهء‌صلمط۴ 
0 ۷ 

۶ )1۳۱۵۵6 :مصیلنناه 18012160 متتعصعع ۵ 10۲ ممتاهاتامه۷-وومی ۵۶ ممتامانصته نان .2012 .م۱001ظ .ظ 24 ریک رلطهمحصه۴ 
34-۰ :53 )صمحصممنا ۲ص 4صه عصلاتباظ .وتماممصنههم لمعمتاماتاممهم 

ومعصمدمو )صمتم و02 01 ]165 ,2014 ,1۲2م۵ک ,۱۷۲ 200 رام ) مقصتا0و۲۱28 ۷۲۰ مصتحاف۴! ,۲۱ .ظ راعصتتا0ظ .ظ بیظ ممصه[م۴ 
فوطت0عع۲۵6ظ عوبامط بمززمرط 2 تفص فلقصصتصه. ها 0۵0ص . تععصهت. فققصظر.. صح . عقفقط . تاامصد هه نان طفناهعط 
ما رمصلزممصاعص لمانه‌تتوخض ۵۶ عمصوعععگصمن مطم1)هصم) 18 

وعقحط 20 ع2عط عصتاام۱۷۵0 .2015 ,۵1 ۷۰ .) 20 رتمک ,۷۲ مصتحامک .ظ مفصتامعع۲۱2 ,۱۷۲ رافحتا0ظ .ظ .ظ .۲ ممصهزم۴ 
25-۰ :136 عصلهممص1مطظ فصمماوومنظ .عمتصمصرول 4تبا تمطمتامانام‌صمی مصتفهه عفبامط تمززمته 2 ۵۶ تعتفصهت 

۱ ۷۰ ,1۷۷۵08 .9 یک ر۷۵۵0 1 .ل مرع1]0ظ ۰( .1 ,۲۱۵۲ ۷۷۰ .و رعظ۲۱۳۷2 .9 ,۲ ملک .1 .۲۱ م۷008 یک .6 رما .ظ بآ ملگ ,960 
مط ما ععبامط عملزمعه 0ماهاتامع۷ پالح‌تنهه ه ۵۶ موه حمتاقاتاصه۷ عمط ۵۶ اصمصصعبممصص .2009 ها ۷۷۰ .۲ 20 بتک 
106-۰ :(1) 104 مطلرهمصنعص فصم)وووماظ .فممتاهایصته تمصمتامانام‌ججم مصتوا ممققعو 6۵10 

تفج ٩۵۵608‏ من فنامتید۷ ما 0مامعزهاناو فصمامتطه عملزمعها رم و1۵ ۴۱۵2۲ .1997 .و۷2 .ل 24 بکاصا .ظ ول مقطمصصن 
665-۰ :(5) 13 متبانهته۸ 1 عصاهعصم من 0عنامجرخ .ععتنهتهمعصه) اصم1طاصنه 

مد موه 2003 .عععلز۵ع0 126ومنمصظ معا ماه ۷- ان تصیطمهتمم‌جصع1 ,2003 .260 ۲ 200 ریک ,120 
۰ مزع 1010ظ مه اههد ۵۶ ساعلهه9 مهعتعمصن رعصتاهع]۱۷ 

٩۷۹۱۵۴ 5۲‏ 2۷12۲ مه 10۲ صمتا۷2112 مه ممتاهاصته متقصصتلی ,۱994 .صقافله ۷۰ مه ررحصقاه0 ۷ .۲ ری مامتان هه2 ۲۷ 
صمص امسانه‌تتوخض وه عمممععتممن 94 ب«ممععخ هه ووع0۵۵۵۲ ۷۷۵۲۱۵ 10۴ ]21 .عصفط عصاوم1 

۰ ۷۵۵)012)10۴8 ]ما مصتازمی ما سمل عنج مضه )0۲مقصه مامتاعدم عصتام۷۵0 .1995 ,معجاه۱۷]2 .۲۱ مضه ری ,۱۷۷۵۲۱۵ 
.231-۰ :(1) 38 موه مط] ۵۶ مصمتامع‌عمه۲ ]1 

۵ 0۲01167 210675 م۷۵۱۵ نهر ,2001 46۲۰صنطو: با مه رتکل۵2۵210۷ ۷۰ به ۷ .۲ ر0صط0مهتا ۸ ریگ ,۷۵202۷ 
724-6۰ :(6) 80 9016868 ناو .فممیتصمی اهامممصممتزاجم طققط ملصا 

کاصفاناامم دباممقهع 10۲ صموره ممتاهاتامه۷ [۱۵62 ۵۶ ممتام0002 .2012 ۷۵۰ .2 204 رما .1 رعصه۷ ٩۲.‏ بتاشنا .ظ .9 بلاط 2 
دص مدرگ )ممجصمما ۱ص کامماومش تمصمتمصنم) 18 12۶ 2012 .ممتاهانصصته نان مصلوده موبامط علژمعه صا ۲۵۱۴۵0۷21 
۰ 21021 1010ظ مه آهعبانم‌تبهد۸ ۵۶ «املهه صهععصصم ر( 2[ 
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